Can supersymmetric loops correct the fermion mass relations in SU(5)? 
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We investigate three different possibilities for improving the fermion mass relations that arise 
in grand unified theories (GUTs). Each scenario relies on supersymmetric loop effects alone, 
without modifying the naive Yukawa unification. First, we consider A-terms that follow the usual 
proportionality condition. In this case SUSY effects can improve the mass relations, but not 
completely. Interestingly, imposing Yukawa coupling unification for two families greatly constrains 
the range of parameters in the MSSM. Secondly, we employ a new ansatz for the tri-linear A-terms 
that satisfies all experimental and vacuum stability bounds, and can successfully modify the mass 
relations. Finally, we investigate the use of general (non-proportional) A-terms, with large off- 
diagonal entries. In this case flavor changing neutral current (FCNC) data present an important 
constraint. We do not pretend to present a complete, motivated theory of fermion masses. Rather 
this paper can be viewed as an existence proof, serving to show that Yukawa coupling unification 
can occur even within the framework of minimal GUTs. 



I. INTRODUCTION 

Unified theories that incorporate supersymmetry have 
met with great success ffl. Included in the attractive 
features of such theories are: the possibility of 
ameliorating the hierarchy problem, a natural radiative 
mechanism for electroweak symmetry breaking, and a 
candidate for dark matter. The most notable success 
of the Supersymmetric Standard Model (MSSM) is the 
successful unification of the gauge couplings . 

In spite of all these attractive features, the minimal 
models appear to be far from complete. Most notably, 
the flavor sector of such theories remains a puzzle. One 
may decrease the number of mysterious parameters in 
the flavor sector by placing quarks and leptons in the 
same multiplets, as required in Grand Unified Theories 
(GUTs). This reduces the number of independent 
Yukawa couplings at the GUT scale However, 
the predictions of such a scheme have met with 
mixed success. For instance, in the minimal SU (5) 
implies the tree-level relation = m/ i at the GUT 
scale, which seems to work quite well for the third 
generation. However the situation for the lighter families 
appears catastrophic, since we would expect the relation 
TO e/ m ^ = Wrf/rTig, but this is off by almost one order of 
magnitude. 

This problem has led some authors to study non- 
minimal models or even to question the whole 
idea of grand unification B. One way to solve this 
problem is to invoke higher-dimensional operators [Q, 
an approach that essentially postpones the problem to 
energies beyond Mqut- Another approach involves the 
use of non-minimal representations for the Higgs sector. 
For instance, in the context of a SU(5) GUT, by enlarging 
the Higgs sector with a 45 representation, Georgi and 
Jalrskog B obtained the relations rrid = 3m e and m s — 



which correctly reproduces the observed fermion 
masses when evolved down to low energies. 

In this paper we present a simpler approach which 
has been overlooked thus far. Instead of adding extra 
matter or higher-dimensional operators to the theory, we 
work within the MSSM framework. We use a radiative 
mechanism to correct the light generation mass relations, 
i.e. — m e , m s — m^, in a manner that respects 

the successful relation m& = m T . It is an exciting 
possibility that one could accommodate the unification 
of the Yukawa couplings within the simplest of SUSY 
theories, the MSSM. 

The correction of fermion masses by SUSY loops has 
been known for a while ||. These corrections have been 
applied in an attempt to generate fully the quark and 
lepton masses H[To). The idea here is to use these 
radiative effects in a somewhat more modest way, to 
fix up the GUT relations, a task which we find to be 
non-trivial. We will investigate this idea within three 
contexts. First, we investigate the case of proportional 
A-terms. Then, we examine two cases with non-minimal 
A-terms. In our second approach, we use a particular 
ansatz that achieves unification for all three families. 
In the final approach, we recognize that A-terms that 
are slightly more generic than the proportional case can 
also help remedy the mass relations. That is to say, 
if one allows off-diagonal A-terms, one can remedy the 
Yukawa relations. It appears that such a scenario is still 
consistent with flavor-changing neutral current (FCNC) 
bounds. 



II. STANDARD A-TERM CORRECTIONS 

We are interested in evaluating the corrections to 
fermion masses coming from SUSY loops. First, we 
consider the standard case where we have A-terms 
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that are proportional to the Yukawa matrices. The 
dominant corrections are the flavor conserving (FC) 
gluino and Higgsino-mediated loops. We include only 
these corrections in our analysis. Since these corrections 
depend on only a few MSSM parameters it is feasible 
to perform a detailed numerical analysis. Because these 
corrections also modify the CKM mixing matrix, we 
prefer to follow a procedure consistent with all low-energy 
data. Our chosen method is essentially a bottom-up 
approach, and consists of these steps: 

• Working in a basis where lepton and u-type quark 
mass matrices are diagonal, we obtain the non- 
diagonal d-quark mass matrices (md)ij, at the 
M^-scale. If we assume that the mass matrix is 
symmetric, we can write: 



m d = Vckm ■ m 



diag 



■v, 



GKMi 



(1) 



where Vckm denotes the quark mixing matrix. 
The lepton mass matrix, mi, is diagonal. We can 
extract the quark mass matrix at Mz utilizing the 
quark masses in Table I, by taking into account the 
QCD renormalization effects. 

• We then modify the mass matrix of the down- 
type quarks by subtracting the contributions of the 
SUSY loops. The SUSY corrections to the lepton 
mass matrix are negligible. The flavor conserving 
(FC) corrections to the down quarks includes the 
gluino and Higgsino loops, namely (Sm d ) FC — 
{6m d y + (Sm d ) 6 . 

• The Yukawa couplings are simply related to the 
mass matrix through the vacuum expectation value 
of the down-type Higgs. Namely, 

m d 



i sin P ' 



(2) 



where v=174 GeV. We convert the Yukawa 
couplings from the MS scheme to the DR scheme, 
using the dictionary of Martin and Vaughn [pTf . 
This is necessary to be compatible with the 
renormalization group equations used in the next 
step. 

Once we have the Yukawa couplings of both quarks 
and leptons at Mz in hand, we evolve the Yukawa 
matrices up to the GUT scale (Mqut), using the 
renormalization group equations of the MSSM [p^[ . 
This yields the Yukawa matrices at the GUT scale 
Y d {M GUT ) and Y l (M GUT ). 

The next step is to diagonalize the Yukawa 
matrix, Y d (McuT), and compare its eigenvalues 
with the lepton Yukawa couplings at the GUT 
scale, Y 1 (Mcut), to see whether unification has 
occurred. To quantify this unification we define 
ei = m-Yj i )/Yi l \. 



Input 


Value (MS scheme) 


a a (M z ) 


0.119 


sin 9w 


0.23124 


<* E m(M z ) 


127.943 


m t (m t ) 


165 GeV 


m b (m b ) 


4.25 GeV 


m c (m c ) 


1.15 GeV 


m s ( 2 GeV) 


115 MeV 


m d ( 2 GeV) 


6 MeV 


m u ( 2 GeV) 


3 MeV 


m T (m T ) 


1.777 GeV 


m M (m M ) 


105 MeV 


m e (m e ) 


0.511 MeV 



TABLE I. Inputs for the numerical analysis of section II. 



• We repeat the above steps, scanning the parameter 
space of the MSSM, to search for those regions 
where q is close to zero. 

Now that we have outlined the steps in the analysis, 
we note the explicit corrections to the Yukawa couplings. 
The gluino loop, (8m d y , contributes to all the diagonal 
elements m ddtBSy bb- We work in the limit when /itan/3 
is much larger than the contribution from the tri-lincar 
scalar coupling. In this limit, the gluino loop contribution 
may be expressed as [Djfl: 



(5m d ) 



a _ 



2a 



flMg 



3tt M\ 



I{x) tan (3m di 



(3) 



The Higgsino-mediated loop (Sm d ) H , is only relevant 
for the third family due to the suppression from small 
Yukawa couplings. Its contribution is given by: 



{Sm d )fj = -jz^Hx) t&n/36 i3 6 j3 . 



(4) 



As usual, Mg denotes the gluino mass, /j, is the parameter 
for the Higgs term in the superpotential, and A t denotes 
the tri-linear scalar term for top. Mx denotes the largest 
mass appearing in the loop, and I(x) denotes the loop 
integral, given by 



I{x) 



-1 + x — x log(x) 

o 3 ^ ■ 



(5) 



where x = rn~/rir~ for the gluino loop, and x = p r/wM 
for the higgsino loop. 

The key point in this analysis is to recognize that the 
higgsino and gluino pieces come with an opposite sign. 
As a result these two contributions may nearly cancel 
each other for the third generation, thereby preserving 
the raj = m T relation to some extent. However, since the 
higgsino piece only contributes to the third generation, 
the gluino piece remains uncancelled for the other two 
families and can be utilized to fix up a mass relation. Of 
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course, it can only be used for one of the two families. It 
will exacerbate the problem for the remaining family. 

The results for the parameter Ci for the 3 generations 
are shown in FIG. 1-3. These graphs correspond to 
m g = mq = 500 GeV, and ,u = -600 GeV. Notice 
that there is a limited region in the (A t , tan 0) plane 
where unification (say |e| < .2%) occurs for both the 
the first and third generation. In fact, if we assume that 
this approach is part of the solution to remedying the 
mass relations, we can view the results of this analysis 
as making predictions for tan (3 and A t . We simply 
mention here that the other sign of /i allows one to 
achieve unification for the second and third generation. 
Though this appears to be more difficult, and requires 
large values of tan/3 where the peturbativity of the b- 
quark Yukawa coupling is an issue. 

One could view this scenario with some wariness, 
as it remedies the mass relation for one generation at 
the expense of exacerbating it for the other generation. 
Nevertheless, this sort of scenario allows one to get the 
correct fermion spectrum by, for instance, introducing 
just one higher-dimensional operator involving the 
remaining family. In some sense, we have achieved "flavor 
from no flavor" . 
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FIG. 1. The percentage deviation in Yukawa couplings, 
100 xe, e 100 x \(Yfi - Yii)/Yji\, for the third generation. 
The dark shaded region represents the area where the b 
quark Yukawa coupling is non-perturbative, and the analysis 
breaks down. The figure corresponds to m g — 500 GeV and 
H = -600 GeV. 



III. A NON-MINIMAL ANSATZ FOR THE 
yl-TERMS. 

Since our initial efforts with proportional A-terms did 
not succeed in achieving complete unification, we are 
forced to resort to a more general form for the A-terms, 



tan (5 

FIG. 2. The percentage deviation in Yukawa couplings, 
100 x €i = 100 x \{Y?i - YjA/Yii], for the second generation. 
The dark shaded region represents the area where the b 
quark Yukawa coupling is non-perturbative, and the analysis 
breaks down. The figure corresponds to m g — 500 GeV and 
M = -600 GeV. 

Percentage Deviation in 1 st Generation Yukawa Couplings 
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FIG. 3. The percentage deviation in Yukawa couplings, 
100 xe, e 100 x \(Yii - Yli)/Yii\, for the first generation. 
The dark shaded region represents the area where the b 
quark Yukawa coupling is non-perturbative, and the analysis 
breaks down. The figure corresponds to m g — 500 GeV and 
ji = -600 GeV. 



thereby opening the door for a larger number of free- 
parameters. This makes it more difficult to perform a 
detailed numerical analysis. Instead, we shall be content 
to present an "existence proof of the viability of these 
mechanisms. In the following, we shall employ a top- 
down approach, by evolving the Yukawa couplings from 
the GUT-scale to the weak-scale, and then including the 
threshold corrections. 

We shall begin by considering a particular ansatz for 
the tri-linear terms, aimed at correcting the Yukawa 
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matrices for the light generations. In this section, we 
say nothing about b-r unification, since that relation 
is already close, and presumably can be fixed up by 
methods similar to those used in the previous section. 
So, let us start by writing the down-type fermion mass 
matrix of first and second generations as: 



M° = 



Am 
Am m 



in 



(6) 



and 



Thus, at GUT scale one has m s 
A 2 = m^/irvl = mP d /rnP s . The superscript will indicate 
fermion masses evaluated at the GUT scale. Now, we 
shall include the correction coming from the gaugino- 
sfermion loops, using the an ansatz for the A-terms that 
results in the following contribution to the squark mass 
matrices: 



l LR 




X 



v cos f3 



A 
A 1 



AY s v cos (7) 



Here Y s is the Yukawa coupling for the strange quark. 
We have separated out the first piece to emphasize the 
largeness of the tri-linear term, X, that is necessary as 
we will see below. The grand unification forces the term 
X to be common to both charged leptons and down- 
type quarks. We note that in the mass-eigenstate basis 
for the fermions, the first matrix on the right-hand side 
becomes off-diagonal, which can lead to FCNC problems. 
We will address this issue shortly. Despite our concerns 
about FCNCs, and the fact that we do not yet have a 
microscopic theory to motivate our ansatz, we will boldly 
forge ahead. We now show that this form for the A- 
terms proves to be phenomenologically very interesting 
in repairing the mass relations. 

The dominant correction to the quark masses arises 
from the gluino-mediated loop. This is larger than the 
bino-loop contribution due to the difference in the gauge 
couplings. Therefore, for questions of unification, we can 
again neglect the corrections to lepton masses. Moreover, 
because of the form chosen for the tri-linear coupling 
terms, Eqn. (|7]), the gluino-loop contributes mostly to the 
(2,2) entry of the quark mass matrix. For simplicity, we 
neglect all corrections but this one. After this correction, 
the lepton mass matrix approximately retains the from 
of Eqn. (o), while the quark mass matrix becomes: 



Am 

Am m 4- Sm s 



(8) 



The basic idea will be to utilize the gluino loop 
correction, Sm s , to cancel a substantial portion of the 
muon mass, m, derived from the Yukawa coupling. 

To determine the precise amount of cancellation 
necessary, one must look at the experimental data at low 
energy. In particular, given the mass matrix of the form 
in Eqn. (JsJ) , one generates a Cabibbo angle, A c given by 
the relationship: 



A 



m d 
to. s 



(9) 



We note that the masses in this equation are the physical 
masses, which yields A c « 0.22. This value is consistent 
with other determinations of A c , a well-known fact. 
Now, utilizing the diagonalized form of Eqn. (||) for the 
corrected mass matrix, we can write: 



A, 



A? 



(mf + Sm s )/(1 + Xl)' 



(10) 



We also note that the ratio A 2 = m^/m^ is independent 
of scale (to one loop), so we can fix the value of A at 
the weak-scale using the physical lepton masses, namely 
A ~ 0.07. This gives us the relation: 



1 



8m s 



(11) 



Also, we note that, after RGE running, we have the 
approximate relation, mf « 2.1m M . 

Finally, we wish to determine the constraint on the 
SUSY parameters. To do this, we need to relate the 
SUSY parameters to Sm s . The correction due to the 



gluino loop can be written as 1 1 3 1 



3tt 



M\ 



(12) 



where m 2 LR is given in Eqn. (0). We assume that m 2 LR 
is completely dominated by the piece proportional to 
X. Here, Mx denotes the largest mass appearing in 
the loop, and I(x) denotes the loop integral, given in 
Eqn. (||), which becomes close to 1 in magnitude at 
small x. Indeed, we will be considering small x, as we 
will see m? ps m? > M|, is necessitated by the FCNC 
data coupled with concerns about naturalness. In the 
following, we denote the mass of all heavy scalar particles 
by m. 

Therefore, combining Eqs. ([?]), (pi]), and (|l2|), we find 
that the supersymmetric parameters must satisfy the 
following constraint if we are remedy the mass relations: 



MgX cos /3 



0.035. 



(13) 



As we intimated earlier, it is a non-trivial task 
to reconcile this requirement with FCNC data. In 
particular, a strong constraint comes from the fi — > 
limit JL4| |. We find that once a gluino mass is chosen, 
Eq. (|13[), along with the limit from /i — > ej, sets the 
mass scale for the sleptons. For the gluino, we choose a 
mass of 600 GeV, and utilize the relation Mg w 1M\ that 
arises from the unification of gaugino masses. A rotation 
of the matrix of Eq. (|j) to the diagonal fermion mass 
basis generates a 8\i in the lepton sector given by 
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$ 2 ) 



LR 



X cos /3 



Xv = 



MgX cos (3\ Xv 
~M~h 



= 7.1 x 10" 



(14) 

Here, we have inserted the value for A and used Eqn. 
(|T3|). Using this value for (S[ 2 )lr and the constraint on 
the [i — ► e7 branching ratio, one can find m ~ 4.4 TcV. 

It then remains to check whether a scenario with 4.4 
TeV sleptons and squarks satisfies the remaining FCNC 
data. Note that the expression for (Sf 2 )LR is identical 
to Eq. (|l4|), with the replacement A — > A c . Using the 
analysis of |L5), which includes NLO QCD effects, we 
have checked that such a scenario satisfies the bounds 
from AMx by about a factor of 5. However, constraints 
from e and e'/e [^6|-^8|, are not satisfied unless the 
phase of (Sf 2 )LR is less than about 0.2. We have also 
checked, adapting the analysis of Jl9|], that this scenario 
is safe from measurements of the anomalous magnetic 
moment of the muon po| . 

Incidentally, considering the large value of the tri- 
linear coupling required in this scenario, one might be 
concerned about charge and color breaking minima of the 
scalar potential. For the values of M~ g and m described 
above, with tan/3 ~ 2, we have X ~ 2500. Although it 
is possible that such a scenario might not have the the 
correct electroweak symmetry breaking (EWSB) vacuum 
as an absolute minimum, X ~ rh does allow the EWSB 
vacuum satisfying the less stringent criteria of meta- 
stability pl| . That is to say, we succeed in producing 
a EWSB vacuum possessing a lifetime greater than the 
age of the universe. 

Finally, we wish to comment on proton decay in this 
scenario. One might expect the proton decay rate to 
be significantly affected on two counts. First of all, 
the first and second generation scalar masses are pushed 
to several TeV. However, increasing the masses for the 
first and second generation superpartners is not enough 
to decrease the proton decay rate predicted by SU(5) 
below experimental bounds, a point implicit in |22] , and 
later emphasized in reference M|. This is due to the 
presence of an important contribution arising from the 
exchange of third generation superpartners that remains 
even if the first two generations decouple. Secondly, 
one might be concerned that the dimension 5 operators, 
which are proportional to Yukawa couplings, might be 
significantly modified. However, since the dimension 5 
operators are proportional to the (2, 1) element of the 
Yukawa matrix, our large correction to the (2, 2) element 
is not significant. Finally, one might be interested 
in entirely new contributions to proton decay arising 
from the new large A-term that might somehow effect 
a cancellation. However, we have checked that such 
contributions are at most sub-dominant. So, the proton 
decay rate in this scenario is comparable to the usual 
supersymmetric SU(5) GUTs, which admittedly have 
difficulty accommodating the experimental bounds. 



Finally, it should be noted that the unification of 
fermion masses with this ansatz requires both scalars of 
a few TeV, while gauginos may remain relatively light. 
Many possibilities exist for such scenarios ^3. However, 
none of these scenarios generates a large enough tri-linear 
coupling, X. Until a micro- physical motivation can be 
found for this case, the above discussion must be simply 
regarded as an existence proof. 



IV. CORRECTIONS FROM 
FLAVOR-CHANGING GLUINO INTERACTIONS 

Finally, let us consider a scenario in which we 
allow effects from flavor violating (FV) gluino-quark- 
squark interactions to enter into our analysis. Such 
contributions are known to induce potentially dangerous 
FCNC effects. Therefore, we must be careful to satisfy 
such bounds. 

In this case, the mass-insertion method proves some- 
what complicated, as several insertions are necessary. 
Therefore, we find it convenient to write the vertex g- 
dh,R-dL,Ri in terms of 3 x 3 mixing matrices, Wl,r — 
V£ R V LR , where V£ R (V L ' a ) denote the rotation matrices 
required to diagonalize the fermion (sfermion LL,RR) 
mass matrices. A discussion of the bounds on the Wl,r 
matrices is presented in JTo| . 

In general, the corrections to fermion masses will 
include corrections from both gluino and Higgsino 
mediated loops. In the large tan/3 limit, it is reasonable 
to include the Higgsino loop only for the heavier fermions 
(third generation), whereas the gluino flavor conserving 
(FC) loop contributes to all diagonal elements of the 
fermion mass matrices. The size of the gluino loop 
is proportional to the fermion mass itself. Thus, we 
can write the correction to quark masses as: Srrid — 
(Smd) FC + {5rrid) FV ■ For the quarks, (Smd) FC is given 
in the previous section in Eqn. (|l2|), the FV gluino 
contribution is given by JIo[| , 

(j, tan/3x 1 / 2 H r 



(5m d ) 



FV 



0.7m;, 



Mg 60 



0.5 



-T. 



3; 



(15) 



where 



L,3i yv R,3i> 



H denotes a combination of 



loop integrals, analogous to Eqn. (^|), and in general H ~ 
0.5. Note that this equation is proportional to mi,, which 
will allow these corrections to reach substantial values, a 
prerequisite for correcting the GUT relations. 

In order to obtain the required values for the fermion 
masses, it seems convenient to fix the values of the 
FC and FV contributions, in such a way that their 
contribution to cancels. Again, we neglect the 
contribution to the leptons. In this case, the ratios 
nidi/m^ can be written as: 

Sm di \ m° di 



= 1 



m di 



(16) 



5 



Here, the superscript again denotes evaluation at 
the GUT-scale. In the above expression, Sm di includes 
contributions from both FC and FV pieces. In order to 
generate the correct mass ratios, 5m di /m di must in turn 
satisfy: 



Sm d \ 
m d J 



FC 



/ X \ p V 

I dm d \ m b 



FC 



( 5m s 



FV 



= 2 



(17) 



(18) 



After substituting the expressions for the factors 
resulting from Eqn. (8), one finds that the correct 
mass ratios require: T^T^ 2 ~ 1CP 5 (for tan/3 = 45). 
This solution is in agreement with current bounds on 
FCNC phenomenology, which imply: T^T^ < 2 x 1CT 5 
]To|J*| Thus, the quark-lepton mass relations for all three 
families can be remedied, using a large effect from both 
FC and FV corrections. 

There could be some concern about whether the 
Yukawa coupling of the b quark would remain pertur- 
bative in such a large tan (3 scenario, or whether this 
solution could survive after more stringent bounds on 
FCNC are applied. In the face of such objections, it 
should be noted that including the FV corrections to the 
lepton masses (which will be proportional to the tau mass 
and could conceivably be non- negligible), or by scanning 
the relevant parameters, one should be able to retain this 
(or a similar) solution. 



while retaining the successful relation mb = m T . This 
ansatz is barely consistent with current FCNC bounds. 

As a final approach to modify the mass relations, 
we considered Flavor violating (FV) gluino-quark-squark 
interactions. The values of the A-matrix required barely 
satisfy all phenomenological and stability bounds of the 
MSSM. 

Although this paper does not present a theory of 
fcrmion masses, it does presents three cases that can be 
considered as "counter-examples" to the usual claim that 
the SU(5) GUT model is excluded because of the wrong 
fermion mass relations. We stress our result: fermion 
mass relations can work for all three families even in the 
minimal SU(5) GUT. 
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V. CONCLUSIONS 

The main result of this paper is that the "wrong 
fcrmion mass relations" of the minimal SUSY GUTs can 
be greatly ameliorated by employing radiative corrections 
due to the superpartners of the SM particles. 

Employing the conventional A-terms, we can obtain 
correct unification relations for the Yukawa couplings for 
two families, depending on the sign of the fi parameter. 
We call this scenario "flavor from no flavor". This 
scenario effectively limits the parameter space of the 
MSSM where two families correctly unify. However, to 
unify the remaining family a higher dimensional operator 
(or a more complicated higgs sector) is required. 

We have also presented a non-minimal ansatz for the 
A-terms, that fixes the relations m d — m e , m s = rrin, 



Tn fact, there are no direct bounds on each factor T3132; 
however, by combining present bounds on several FCNC 
processes involving the R elements, one can arrive at the 
quoted bound 
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